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Abstract 

Quasi-wavelets  (QWs)  are  a  representation  of  turbulence  consisting  of  self-similar,  eddy¬ 
like  structures  with  random  orientations  and  positions  in  space.  They  are  used  in  this  paper  to 
calculate  the  scattering,  due  to  turbulent  velocity  fluctuations,  of  sound  behind  noise  barriers 
as  a  function  of  the  size  and  spatial  location  of  the  eddies.  The  sound  scattering  cross-section 
for  QWs  of  an  individual  size  class  (eddy  size)  is  derived  and  shown  to  reproduce  results  for 
the  von  Karman  spectrum  when  the  scattered  energies  from  a  continuous  distribution  of  QW 
sizes  are  combined.  A  Bragg  resonance  condition  is  derived  for  the  eddy  size  that  scatters  most 
strongly  for  a  given  acoustic  wavenumber  and  scattering  angle.  Results  for  scattering  over 
barriers  show  that,  for  typical  barrier  conditions,  most  of  the  scattered  energy  originates  from 
eddies  in  the  size  range  of  approximately  one-half  to  twice  the  size  of  the  eddies  responsible  for 
maximum  scattering.  The  results  also  suggest  that  scattering  over  the  barrier  due  to  eddies 
with  a  line  of  sight  to  both  the  source  and  receiver  is  generally  significant  only  for  frequencies 
above  several  kilohertz,  for  sources  and  receivers  no  more  than  a  few  meters  below  the  top  of 
the  barrier,  and  for  very  turbulent  atmospheric  conditions. 
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1.  Introduction 

Turbulence  can  reduce  the  effectiveness  of  outdoor  noise  barriers  by  scattering  sound 
energy  deep  within  the  diffractive  shadow  region  (that  is,  the  indirectly  insonified  region) 
formed  by  the  barrier.  Accurate  modeling  of  this  phenomenon  is  quite  challenging,  as 
realistic  representations  for  the  turbulence  in  the  vicinity  of  the  barrier  are  required  in 
addition  to  techniques  for  calculating  the  complex  sound  propagation  effects.  Daigle  [1] 
provided  the  first  detailed  calculations  of  the  scattered  and  diffracted  sound  energy 
behind  a  barrier,  and  compared  the  results  to  experiments  involving  detailed  acoustical 
and  turbulence  measurements.  He  demonstrated  that  scattering  can  significantly  impact 
barrier  performance,  particularly  for  frequencies  above  a  few  kilohertz. 

Recently,  a  new  tool  for  representing  turbulence,  called  quasi-wavelets  (QWs) 
[2,3],  was  introduced.  Unlike  Fourier  spectral  representations  of  turbulence,  QWs 
have  the  advantage  of  being  localized  in  space.  This  feature  makes  scattering  cal¬ 
culations  more  straightforward  and  allows  the  scattering  process  to  be  characterized 
simultaneously  in  space  and  scale  of  the  scatterers.  Another  improvement  is  that  the 
QW  representation  includes  scattering  from  the  energy-containing  subrange  of  the 
turbulence  (the  large  eddies).  Goedecke  et  al.  [4]  successfully  used  a  QW  formulation 
to  calculate  the  spectral  broadening  of  sound  scattered  by  moving  eddies.  In  this 
paper,  we  reexamine  the  phenomenon  of  sound  scattering  over  barriers  using  a  QW 
representation  of  the  turbulence. 

The  individual  QWs  are  intended  to  have  properties  reminiscient  of  actual  eddies 
in  a  turbulent  flow.  The  general  technique  of  representing  turbulence  in  wave  scat¬ 
tering  calculations  as  a  collection  of  discrete,  eddy-like  structures  (sometimes  called 
turbules)  apparently  originates  with  DeWolfc  [5].  More  recently,  McBride  et  al.  [6] 
used  such  a  representation  for  scattering  of  sound  above  a  complex  impedance 
boundary,  and  Boulanger  et  al.  [7]  used  it  to  determine  turbulence  effects  on  sonic 
booms.  We  call  the  eddy-like  structures  quasi-wavelets  when  the  ensemble  is  con¬ 
structed  from  a  parent  function  in  a  systematic,  self-similar  manner. 

The  paper  is  organized  as  follows.  We  provide  an  overview  of  the  QW  model 
formulation,  together  with  some  new  results  for  determination  of  the  model  pa¬ 
rameters,  in  Section  2.  Equations  are  derived  in  Section  3  for  the  scattered  field  from 
a  single  QW  and  for  the  scattering  cross-section  of  a  random  spatial  distribution  of 
QWs  belonging  to  a  single  size  class  (eddy  size).  When  the  scattered  energies  from  a 
continuous  distribution  of  QW  size  classes  are  combined,  the  resulting  cross-section 
agrees  with  a  previous  result  for  the  von  Karman  turbulence  spectrum.  In  Section  4, 
QWs  are  applied  to  the  problem  of  sound  scattering  over  a  barrier.  The  results  il¬ 
lustrate  the  eddy  scales  primarily  responsible  for  scattering  as  well  as  the  spatial 
regions  from  which  the  scattered  energy  originates. 


2.  Quasi-wavelet  ensembles 

Turbulence  is  generally  conceived  as  a  collection  of  eddies  of  many  different  sizes 
[8,9].  The  QW  model  is  a  mathematical  representation  of  turbulence  that  more 
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closely  resembles  this  physical  picture  than  Fourier  modes  or  customary  wavelets. 
Like  customary  wavelets  [10,11],  the  QW  representation  is  based  on  self-similar 
application  of  a  spatially  localized  parent  function.  In  contrast  to  customary 
wavelets,  however,  the  orientations  and  positions  of  the  quasi-wavelets  are  random. 
Furthermore,  the  QW  basis  functions  are  not  required  to  be  mutually  orthogonal  or 
to  form  a  mathematically  complete  set. 

Recently,  a  QW  representation  has  been  developed  that  reproduces  the  von 
Karman  turbulent  velocity  spectrum  [3],  which  generally  fits  atmospheric  turbulence 
well.  In  [2],  expressions  for  temperature  and  velocity  fluctuation  spectra  for  isotropic 
homogeneous  turbulence  were  obtained  using  QW  (turbule)  superpositions.  Only  the 
velocity  spectra  are  considered  in  this  paper.  The  relationship  between  the  QW 
parent  function  and  the  turbulence  kinetic  energy  spectrum  is  discussed  in  this  sec¬ 
tion  in  detail.  The  parameters  of  the  QW  ensemble  are  related  to  the  inertial- 
subrange  characteristics  of  the  turbulence. 

2.1.  General  results 

Like  customary  wavelets,  QWs  are  derived  from  translations  and  dilatations  of  a 
dimensionless,  spatially  localized  parent  function  /{Cj.  Here,  ^  is  the  magnitude  of 
the  vector  ^  =  (r  -  b“”)/a„,  where  r  is  the  spatial  coordinate,  b“”  is  the  center  of  the 
QW,  and  is  its  size.  The  index  a  indicates  the  size  class  of  the  QW,  with  a  —  1 
being  the  largest  size  and  a  — N  the  smallest,  and  «  indicates  a  particular  QW  within 
that  size  class.  The  size  ax  is  associated  with  the  outer  scale  of  the  turbulence  and  a^] 
with  the  inner  scale.  By  defining  the  velocity  field  v“"(r)  as  the  curl  of  a  vector  po¬ 
tential  A“”(r),  we  assure  that  the  turbulence  will  be  a  solenoidal  field.  Hence,  we  set 
v“"(r)  =  V  X  A“"(r),  where 

A“"(r)  =  n“''a^(a  (1) 

and  is  the  angular  velocity  vector  of  the  QW.  The  presence  of  a^  in  the  definition 
provides  dimensional  consistency.  Writing  out  the  curl  of  the  potential  leads  to  the 
following  result  for  the  rotational  velocity  field  associated  with  the  QW: 

v""(r)  =  X  (r  -  b““)(  -  r'e//e^).  (2) 

Since  turbulence  and  scattering  processes  are  often  studied  in  spectral  domain,  we 
will  have  use  for  the  Fourier  transform  of  the  velocity  field.  By  definition, 

(3) 

{In)  J 

in  which  k  is  the  wavenumber  vector  for  the  turbulence.  Transformation  of  Eq.  (2) 
leads  to 


v“"(k:)  =  iic  X  A“"(k)  =  i(K  x  exp  ( --  Ik  •  \f^)alF{Kaa), 


(4) 
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where 

=  j  (5) 

is  the  spectrally  transformed  parent  function. 

A  key  distinction  between  QWs  and  customary  wavelets  is  that  the  orientation  of 
the  QWs,  ft"*  (where  the  carat  indicates  a  unit  vector),  and  their  eddy  centers,  b“", 
are  random  variables.  For  a  homogeneous,  isotropic  model,  the  Q*"  are  statistically 
independent  and  distributed  uniformly  over  all  directions,  and  the  b®"  are  statistically 
independent  and  distributed  uniformly  in  a  volume  of  interest  V.  This  construction  is 
illustrated  in  Fig.  1.  The  total  velocity  field  created  by  the  quasi-wavelets  (the  QW 
ensemble),  found  by  superposition,  is 

’'('■)  =  (6) 

a=! 

where  N  is  the  number  of  size  classes  and  is  the  number  of  QWs  for  the  size  class 
a.  We  may  define  the  spectral  density  tensor  <^,j(ic)  of  the  velocity  fluctuations  as 

^oW=^(«i.WD;(K)),  (7) 


Fig.  I.  Distribution  of  quasi-wavelets  in  a  volume  V.  Each  QW  has  a  random  position  b"  and  angular 
velocity  vector  ft*",  where  a  is  an  index  for  the  size  class  of  the  QW  and  n  is  an  index  for  the  QWs  within 
each  size  class.  The  nominal  size  of  the  QW  is  a,.  Also  shown  is  the  geometry  used  for  calculating  scat¬ 
tering  by  a  QW  in  Section  3.  The  source  position  is  indicated  as  S,  the  receiver  position  as  R,  and  the  origin 
as  O.  The  vector  r  extends  from  the  origin  to  the  field  point  at  which  the  scattering  is  calculated. 
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where  ( )  is  the  ensemble  average  (average  of  a  large  number  of  random  realizations) 
over  the  H  and  b  variables.  Substituting  with  Eq.  (4)  and  recalling  that  the  QWs  have 
statistically  independent  orientations,  we  find 

<P,(k)  =  M  f^N^a'^F\Ka^){[{K  x  H'")  •  e,]  [(k  x  £1“)  •  e,]),  (8) 

^  a=l 

where  e,  is  the  unit  vector  in  the  ith  direction.  Since  the  angular  velocity  magnitude  is 
the  same  for  all  members  of  a  size  class,  we  may  define  =  |n“"|.  Writing  out  the 
components  of  this  equation  and  noting  that  the  cross-correlations  of  the  angular 
velocity  components  (e.g.,  are  zero,  whereas  the  autocorrelations  (e.g., 

((0“”)^))  must  equal  0^/3,  we  find 

<?y(K)  =  -  KiKj)  53 N.i^yj’F^Ka,),  (9) 

a=l 

where  dij  —  1  when  i  =  j  and  0  otherwise.  As  is  well  known  in  the  theory  of  isotropic 
turbulence,  the  spectral  tensor  can  be  written  as  =  E{K){K^dij  -  KiKj)j{AnK^), 
where  E{k)  is  the  turbulent  kinetic  energy  (per  unit  mass)  spectrum  [9].  Using  this 
formula  and  Eq.  (9),  we  have 

E{K)  =  ^^Y.^.ny:F\Ka,).  (10) 

a=l 

The  spectral  properties  of  the  QW  model  are  dependent  upon  the  choice  of  F(y)  and 
the  scaling  of  and  with  Let  us  construct  the  representation  such  that  the 
packing  fractions, 

ci>,=N,c^JV,  (11) 

equal  a  constant  value  cp,  independent  of  the  eddy  size  class  a.  Next,  we  choose  the 
angular  velocity  magnitudes  according  to  Kolmogorov’s  [12]  original  hypotheses  for 
the  inertial  subrange.  According  to  that  treatment,  quantities  in  the  inertial  subrange 
depend  only  on  the  eddy  size  (a^  in  our  notation)  and  the  dissipation  rate  of  specific 
turbulent  kinetic  energy,  €.  By  dimensional  analysis  (since  e  has  dimensions  length^ 
time“^),  we  must  have 

=  (12) 

where  cq  is  a  constant.  Substituting  these  scalings  into  (10)  yields 

(13) 

^  a=l 

This  result  shows  that  we  can  adjust  the  product  (pc^  in  order  to  obtain  a  desired 
turbulent  kinetic  energy. 
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As  with  customary  wavelets,  the  QW  size  classes  are  chosen  to  scale  in  a  self¬ 
similar  fashion;  that  is,  —  const.  For  convenience,  we  define  a  constant  //  >  0 

such  that 


(14) 


When  the  size  classes  are  very  closely  spaced,  this  relationship  implies  = 
«a+i  -  fla  —  Therefore, 


E{k)  ~ - 

3, 


With  the  substitution  y  —  kq,  we  have 

32T:>4e2/3K-5/3  r-.. 

— L  ^ 


E(k) 


(15) 


(16) 


The  ratio  (p/^i  can  be  removed  from  this  expression  by  noting  that  the  total  number 
of  quasi-wavelets  in  the  representation  is 


N  N 

Aqw  =  ^ 

a=l  «=I  daw 


■\  a=l 

The  integration  results  in 

^ow  ^  K’  -  07’) 


Solving  for  the  ratio  cp/fi,  we  may  now  rewrite  Eq.  (16)  as 
327t^A^Qw4€'''5K-V5  ,  r-i 

—  (V-"i  )  , 

xaw 


E{k)  ~ 


(V  -"1’) '  /  d>>y'‘/^F-(v). 


(17) 


(18) 


(19) 


To  maintain  a  fixed  value  of  the  energy  spectrum  E{k),  must  be  varied  in 
proportion  to  the  inverse  of  the  number  of  quasi-wavelets  per  unit  volume. 

We  next  turn  to  the  matter  of  relating  the  parameters  of  the  QW  ensemble  to 
known  properties  of  the  turbulence.  The  energy  spectrum  in  the  inertial  subrange  has 
been  shown  empirically  to  obey  the  equation  [13,14] 

£(K')-(55a,yi8)e2/3K-V3^  (20) 

where  a,„  an  empirical  constant,  has  the  approximate  value  0.52.  Comparing  Eqs. 
(19)  and  (20),  we  have  for  wavenumbers  k  the  inertial  subrange  (af'  <  k:  <  a"*), 

5^  =  32^^^^  („-3  _  n 

dxON 

The  preceding  equality  can  be  satisfied  by  setting 


(21) 
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2  5Scc^V  /  _3\  _  5S(Xf)H 

""  5767iWqw  ^  “  192iiV 

and  normalizing  the  function  F(y)  such  that 


pKOl 

/  d>’ 

J  KflV 


Since  kgi  >  1  and  kgi^  <  1  for  k:  in  the  inertial  subrange  (and  anticipating  that 
realistic  spectral  parent  functions  will  be  monotonically  decreasing  in  y,  well  behaved 
near  y  =  0,  and  very  small  for  y  >  1),  the  limits  on  the  integral  may  be  replaced  by  0 
and  00,  resulting  in  the  following  normalization  condition: 


roo 

/  dy  y'^'^F'^iy) 
Jo 


Combining  (22)  and  (19),  we  have 


£(^)  =  5^,2/3^-5/3  r 
lo  Jkom 


Eq.  (25)  provides  the  energy  spectrum  from  a  prescribed  F{y).  In  some  applications, 
we  may  wish  to  reverse  the  process  by  calculating  FO)  from  a  prescribed  energy 
spectrum.  Multiplying  both  sides  of  Eq.  (25)  by  putting  =0  (neglecting  the 
effect  of  dissipation  on  the  energy  spectrum),  and  differentiating  with  respect  to  k 
yields  an  equation  for  F^(K:ai), 


F^{Kai)  = 


-2/3^-17/3  -14/3 


[k^^^E{k)]. 


2.2.  Von  Karman  quasi-wavelet 

The  von  Karman  energy  spectrum  has  been  found  useful  for  calculations  of  sound 
propagation  in  the  turbulent  atmosphere  [15,16].  It  is  given  by 


E^{k)  = 


55r(5/6)  oyil 
9i/iir(l/3)  (l  +  fc^IJ)' 


where  is  the  variance  and  Ly  an  outer  length  scale.  Taking  the  limit  KLy  oo  and 
matching  the  result  with  Eq.  (20)  for  the  inertial  subrange  yields  a  relation  among 
Ly^Oyj  and  e: 

_  2r(5/6)  (28) 

[a,y/ir{\/3)\  £' 

A  spectral  parent  function  that  yields  the  von  Karman  spectrum  can  be  deter¬ 
mined  from  Eq.  (26).  The  result  is 
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It  can  be  shown,  by  substituting  Eq.  (29)  into  Eq.  (24),  that  the  normalization 
condition  is  satisfied  for  any  value  of  the  ratio  Lv/oi.  Normally,  one  would  take  the 
largest  size  class  to  be  near  the  outer  scale;  that  is,  L^/a\  ~  1. 

The  spatial  parent  function  f^{^)  (the  “von  Karman  QW”)  is  found  by  substi¬ 
tuting  Eq.  (29)  into  the  inverse  of  Eq.  (5).  Because  of  the  spherical  symmetry,  the 
inverse  transform  reduces  to  the  following  integral: 

m  =J  JJ  dy  ysm{yi)Fiy).  (30) 

With  F{y)  given  by  Eq.  (29),  the  integral  can  be  found  in  tables  [17].  After  some 
algebra,  the  result  is 


/v(0  - 


j  71/227/12 
3J/2r(23/12) 


(31) 


where  Ky  is  the  modified  Bessel  function  of  the  second  kind. 

Fig,  2  shows  QW  energy  spectra  based  on  Eq.  (29)  for  representations  with  =  4 
and  Na  =  64  size  classes.  The  actual  von  Karman  energy  spectrum  (Eqs.  (27)  and 
(28))  and  the  Kolmogorov  spectrum  (Eq.  (20))  arc  shown  for  comparison.  The  length 
scales  for  the  calculations  are  such  that  a\  =  -  lOOOa/v.  We  see  that  the  QW 


Fig.  2.  Energy  spectra  resulting  from  QW  representations  with  4  and  64  discrete  size  classes.  The  von 
Karman  and  Kolmogorov  energy  spectra  are  shown  for  comparison. 
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representation  with  64  size  classes  is  essentially  the  same  as  the  von  Karman  spec¬ 
trum  within  the  energy-containing  {kL^  <  5)  and  inertial  (5  <  kL^  <  500)  subranges. 
Within  the  dissipative  subrange  {kU  ^  500),  the  QW  representation  exhibits  rapid 
diminishment  in  the  spectrum  with  increasing  wavenumber.  This  spectral  roll-ofif, 
which  is  similar  to  that  found  in  actual  turbulence,  is  not  present  in  the  von  Karman 
spectrum  because  Eq.  (27)  does  not  model  it.  The  QW  representation  with  four  size 
classes  produces  a  very  bumpy  spectrum  that  is  not  consistently  close  to  the 
von  Karman  spectrum.  These  results  suggest  that  QW  representations  with  ap¬ 
proximately  20  or  more  size  classes  per  decade  provide  smooth,  accurate  spectral 
models. 


3.  Scattering  by  quasi-wavelets 


Let  us  first  consider  the  scattering  from  a  single  QW.  Fig.  1  defines  the  scattering 
geometry.  Here,  R  is  the  vector  from  the  origin  to  the  receiver,  Rq  is  the  vector  from 
the  origin  to  the  source,  n  =  R/|R|,  no  =  ~Ro/|Ro|»  k  =  An,  ko  =  /mo,  and  6  is  the 
angle  between  n  and  no.  The  acoustic  wavenumber  k  is  27i//co,  where  /  is  frequency 
and  Co  is  the  background  sound  speed.  The  pressure  field  associated  with  the  wave 
propagating  spherically  from  the  source  is  po{R  -  Ro)  =  A  exp(iA:|R  -  Ro|)/|R  -  Ro|, 
where  A  is  the  pressure  amplitude  1  m  from  the  source.  For  single  scattering  by 
velocity  fluctuations,  the  scattered  pressure  field  is  [18,  Eq.  (7.6)] 


2  COS  9  no  ♦  v(r) 


Co 


(32) 


The  primed  vectors  Rq  and  R'^  extend  from  a  point  in  the  scattering  volume  to  the 
source  and  receiver,  respectively.  With  the  approximations  kR[  ~  —  r  •  k  and 

kR^  ~  A:/?o  +  r  •  ko  (valid  for  R,Ro  >  where  Lq  is  the  outer  length  scale 

of  turbulence,  as  discussed  in  section  7,1.2  of  [18]),  we  may  interpret  the  integral  in 
the  previous  expression  as  a  Fourier  transform.  Applying  Eq.  (3),  we  thus  have 

2cos^  no  •  y(-K)  ■ 

where  K  =  ko  -  k  =  /r(no  -  n)  is  the  scattering  vector.  For  a  single,  rotating  QW, 
substitution  with  Eq.  (4)  gives  a  formula  for  the  single-scattered  sound  field 


A  Ir^ 

;%(R)  =  — ^exp[ii(/?  +  /?o)] 


pr(R)  = 


lAn^Ak'^alFiKa^)  cose 
coRRo 


exp  [ik{R  +  Ro-{n-  no)  •  b"’)][no  ■  (K  x  £1”")]. 


(34) 


The  ranges  of  applicability  of  the  single  scattering  approximation  are  given  by  Eqs. 
(7.13)  and  (7.14)  in  [18]. 

The  total  scattered  field  is  a  sum  of  Ps”(^)  ^ 

fields  scattered  by  all  QWs.  To  calculate  the  mean  intensity  of  the  total  scattered 
field,  one  needs  to  multiply  the  total  field  by  its  complex  conjugate  and  then  average 
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over  random  positions  b“”  and  angular  velocities  H*”.  Since  the  b“"  are  statistically 
independent  and  randomly  distributed,  the  contribution  to  the  intensity  will  be  zero 
for  the  cross-terms  involving  two  distinct  QWs  (different  value  of  a  or  n).  Therefore, 
the  intensity  of  the  total  scattered  field  is  the  incoherent  sum  of  intensities  of  the 
sound  fields  scattered  by  the  individual  QWs. 

Let  us  now  determine  the  mean-square  sound  pressure  {|/?s”|^)  for  an  individual 
QW  with  varying  randomly  in  orientation.  This  requires  an  evaluation  of 
([no  *  (K  X  Writing  out  the  components  of  this  equation  and  applying  the 

properties  of  the  angular  velocity  correlations  that  were  described  in  Section  2.1,  we 
find 

((no-(KxJl”)f)=^|n„xiip.  (35) 

Since  |no  x  np  =  sin^ 0 and  AT  =  ky/l{\  -  n-no)  =  k-^2{\  -cos0)  =  2/fcsin(0/2),  we 
thus  have  for  the  mean-square  pressure  field  scattered  by  a  single  QW, 

_  \(>Ti''A^k‘’Qla'^F'^{2ka,  sin  0/2)  sin^  0cos^  0 
/ - 

For  the  von  Karman  spectral  function,  Eq.  (29),  this  becomes 


212iiWQla[^sm^ 


Ocos^  9  f  Lv 


W/J 


(37) 


To  find  the  scattered  mean-square  pressure  from  all  QWs  of  a  given  class  size 
within  some  scattering  volume  K,  the  contributions  from  each  QW  are  simply  added 
incoherently,  as  described  earlier:  {|^|p)  =  E!I'ii(l/’s"P>  Making  sub- 

stitutions  with  Eqs.  (11),  (12)  and  (22),  we  therefore  have 


935oL^,f.iA^k^e^^^  sin^  0cos^  0 

mclR% 


X 


4k^Llalsm\e/2)y^^''' 

) 


(38) 


This  result  provides  a  decomposition  of  the  scattering  based  on  the  size  of  the  tur¬ 
bulent  eddies.  It  applies  to  scattering  from  both  the  energy-containing  and  inertial 
subranges  of  the  turbulence.  In  the  following  section  of  this  paper,  we  will  use 
Eq.  (38)  to  determine  the  location  and  size  of  the  eddies  primarily  responsible  for 
scattering  of  sound  energy  over  barriers. 

The  scattering  cross-section  per  unit  volume  into  a  unit  solid  angle  is  defined  as 
(e.g.,  [18],  p.  188) 
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where  {/s)  =  ([psf)/(2po<^o)  is  the  scattered  intensity  and  /q  —A^/{2pqCqRI)  is  the 
unscattered  intensity  incident  on  the  volume.  Eq.  (38)  leads  to  the  following  scat¬ 
tering  cross-section  for  each  size  class: 


o^{e)  = 


sin^  0cos^  0 


n 


Consider  the  ratio  of  the  cross-section  for  scatterers  of  size  class  a  to  size  class  1. 
Using  the  previous  equation,  we  can  write  this  as 


a' {6)  V"i/  }  +  {qaJaif  _ 


(41) 


where  q  =  2kL,  sin  ( 0/2)  =  kL^  is  a  normalized  scattering  vector.  Fig.  3  shows  this 
ratio  as  a  function  of  q  for  various  values  of  a\/aa.  Note  that  relatively  large  eddies 
scatter  most  efficiently  for  small  q  (forward  scattering  and/or  wavelengths  small 
compared  to  Ly),  whereas  relatively  small  eddies  scatter  most  efficiently  for  large  q 
(large-angle  scattering  and/or  wavelengths  small  compared  to  Ly),  as  previously 


2  k  sin(e/2) 

Fig.  3.  Ratio  of  the  scattering  cross-section  (f  for  QWs  in  the  size  class  a  to  c’  for  QWs  in  the  largest  size 
class.  Results  for  various  values  of  a\/aa  are  shown  as  a  function  of  2A:Lv  sin(0/2),  where  Ly  is  the  von 
Karman  outer  length  scale  and  6  the  scattering  angle. 
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emphasized  in  [4],  By  differentiating  the  preceding  expression  with  respect  to  the 
ratio  ajau  one  finds  that  QWs  with  the  size  ajai  =  y/\l / 6/ q  are  the  strongest 
scatterers.  Alternatively,  we  could  write  this  condition  as  ka^sm{6/2)  = 
yf\l jlA{a\jLy),  which  casts  the  result  as  a  Bragg  scattering  resonance  for  a  given 
value  of  the  ratio  fli/Lv 

To  check  Eq.  (38),  let  us  compare  it  now  to  previous  results  for  scattering  from  a 
continuous  distribution  of  eddy  sizes.  Summing  Eq.  (38)  over  the  size  index  a  and 
then  taking  the  continuous  limit  (as  was  done  in  a  similar  context  in  Section  2.1),  we 
have 


93542^‘/^a,f2/3|Zsin^0cos2  0 


Jka, 


mclR^Rl 

2 


14/3 


I  Vz.;sin  (0/2)  \ 


■23/6 


(42) 


As  before,  we  may  replace  the  lower  limit  on  the  integration  by  0  when  k  is  within  the 
energy-containing  or  inertial  subrange.  The  integration  may  then  be  performed  by 
making  the  substitution  /  =  1  +  (2Lv  sin{^/2)/ai)”^y"^,  with  result 


Fig.  4.  Total  scattering  cross-section  (normalized  by  multiplication  with  cJ^yVC;)  for  the  von  Karman 
turbulence  spectrum  as  a  function  of  scattering  angle.  Various  values  of  JcLy  are  shown. 
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(N-) 


V  sin^  0  cos^  6 

F  ;  ,-|  17/6  • 

36clR^Rl  (2  sin  {O/Tjf  +  {kL,y^ 


(43) 


Substituting  the  structure-function  parameter  Cj  =  (3/2)r(l/3)aye^/^  [16,  Eq.(20)], 
the  overall  scattering  cross-section  becomes 


a{e)  = 


sin^  9  cos^  9 


54r{l/3)cl  [(2  sin  (0/2))'  +  (tiv)'"] 


17/6  * 


(44) 


This  result  agrees  with  [19],  which  was  derived  for  the  von  Karman  velocity  spec¬ 
trum.  For  the  inertial  subrange,  A:Lv  >  1,  the  cross-section  reduces  to  Eq.  (7.17)  in 
[18],  which  was  derived  from  the  Kolmogorov  spectrum.  Fig.  4  shows  the  depen¬ 
dence  of  Eq.  (44)  on  9  for  various  values  of  kLy^. 


4.  Application  of  QWs  to  scattering  over  barriers 

The  geometry  of  the  barrier  scattering  problem  is  illustrated  by  Fig.  5.  For  sim¬ 
plicity,  in  this  paper  the  ground  is  assumed  rigid  and  the  reflections  are,  therefore, 
represented  by  an  image  source  and  an  image  receiver.  We  consider  four  possible 
transmission  paths  involving  scattering  by  turbulence:  source-scatterer-receiver  (line 
SOR),  source-ground  reflection-scatterer-receiver  (line  WOR),  source-scatterer- 
ground  reflection-receiver  (line  SOU),  and  source-ground  reflection-scatterer- 
ground  reflection-receiver  (line  WOU). 


Fig.  5.  Geometry  for  sound  scattering  over  a  barrier.  A  scattering  volume  V  populated  with  random  QWs 
is  shown  above  the  barrier.  S  is  the  position  of  the  actual  source,  W  the  image  source,  R  the  actual  receiver, 
and  U  the  image  receiver. 
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To  examine  the  spatial  properties  of  the  scattering  process,  we  partition  the  total 
scattering  volume  V  into  a  number  of  subvolumes  dVi  and  calculate  the  scattering 
cross-section  o’- (0,)  for  each  subvolume  and  QW  size  class.  As  a  practical  matter, 
each  subvolume  is  treated  as  an  effective  point  scattercr.  This  is  reasonable  when  the 
distances  between  the  subvolume  and  the  source  and  receiver  {R  and  /?o,  respectively) 
are  many  wavelengths  long.  Note  that  the  derivation  of  the  scattering  cross-section 
in  the  previous  section  already  assumed  this  condition.  Neglecting  for  a  moment  the 
ground  reflections,  Eqs.  (38)  and  (40)  allow  us  to  write  the  received  mean-square 
scattered  pressure  attributable  to  the  given  scattering  subvolumc  and  size  class  in  the 
following  form: 

(IpH')  =^'|G(Ro./  +  R*,nRw)l^^f^-T‘'(0,)|G(Rw,R/  -f  R*,,)|', 

where  Rq.,,  R/,  and  R/, ,  are  the  vectors  from  the  center  of  the  scattering  subvolume  to 
the  source,  receiver,  and  the  top  of  the  barrier,  respectively.  The  angle  0,  is  the 
scattering  angle  for  the  subvolume,  as  defined  in  Fig.  5.  The  function  G(Rs,  Rr)  is  the 
Green  function  for  propagation  from  the  source  at  position  Rs  to  the  receiver  at 
position  Rr  with  the  vectors  defined  relative  to  the  top  of  the  barrier.  In  free  space, 
we  would  have  simply  G(Rs,Rr)  =  exp(iA'|R,s  -  Rr|)/|Rs  -  Rr|.  However,  it  is  more 
suitable  to  use  the  Green  function  for  the  sound  field  in  the  vicinity  of  a  barrier  that 
includes  diffraction  behind  the  barrier.  The  calculations  in  this  paper  use  Pierce’s 
[20]  Eq.  (9-9.1),  for  this  purpose.  (The  function  G(Rs,Rr)  here  is  p/S  in  Pierce’s 
notation.) 

Extending  the  free-space  result  to  include  contributions  from  the  reflected  paths, 
we  have 

(brP)  =  [|C(Ro,  +  R*,,  R*,,)|'^f^<T”(0,)|G(R*,-,  R/  +  R/,,,)P 

+  |G{R„,,.  +  R,,„  R,.,)|^5f^a"'(0,  +  ft)|G(R*,„  R,  +  R*,,)P 
+  |G(Ro,,-  +  Rw,  +  x,)|G(R*.,,  R„,,.  +  R*,,)|' 

+  |G(R„,,,-  +  R,,,.,  +  ft  +  j;,)|G(Rw,  R„,  +  R,.,)|']  (45) 

where  R„.,  is  the  vector  from  the  center  of  the  scattering  subvolume  to  the  image 
source,  R„ ,  the  corresponding  vector  to  the  image  receiver,  and  the  angles  ft  and  Xi 
are  defined  in  Fig.  5.  The  four  terms  in  the  equation  correspond  to  the  propagation 
paths  SOR,  WOR,  SOU,  and  WOU,  as  described  previously.  The  scattered  contri¬ 
butions  from  each  of  the  four  paths  are  assumed  to  be  uncorrelated  in  this  formu¬ 
lation.  The  reader  is  referred  to  [21]  for  a  discussion  of  path-correlation  effects. 

Except  as  otherwise  noted,  the  calculations  to  follow  in  this  section  are  based  on  a 
barrier  of  height  2.5  m,  with  Zs  =  0.5  m,  Xs  =  -10  m,  =  1.2  m,  and  jCr  =  15  m.  The 
total  volume  used  for  the  scattering  calculations  spanned  -20  m^x:^25  m, 
-20  m^y<20  m,  and  0  m^z<40  m.  This  overall  volume  was  partitioned  into 
60  X  21  X  60  subvolumes.  In  the  QW  representation,  the  scale  of  the  largest  size 
class,  fli,  is  set  to  5  m.  This  value,  which  is  representative  of  the  size  of  the  largest 
shear-generated  eddies  in  the  vicinity  of  the  barrier  top,  is  also  used  for  Ly.  (In  [16] ,  it 
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Cumulative  contribution  to  scattering  (%) 

Fig.  6.  Vertical  cross-section  through  the  region  of  primary  scattering,  as  determined  from  the  QW  model. 
The  calculation  is  for  a  frequency  of  2000  Hz,  barrier  height  2.5  m,  —  0.5  m,  Xs  =  -10  m,  Zr  =  1  -2  m,  and 
Xr  —  15  m.  Ten  percent  of  the  scattering  occurs  outside  the  10%  contour,  20%  of  the  scattering  occurs 
outside  the  20%  contour,  etc. 


Fig.  7.  Comparison  of  scattering  calculations  performed  with  the  Kolmogorov,  von  Karman,  and  QW 
models.  A  total  of  128  QW  size  classes,  ranging  from  oi  =  5  m  to  =  0.005  m,  were  used. 
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was  shown  that  =  L8z  for  shear-generated  turbulence  over  flat  ground,  implying 
that  Iv  =  5  m  when  z  =  2.8  m.)  The  scale  of  the  smallest  size  class,  is  set  to  0.005 
m,  which  is  characteristic  of  the  Kolmogorov  microscale  in  the  atmosphere.  A  total 
of  128  size  classes  are  used,  in  order  to  obtain  a  very  smooth  spectrum. 

For  simplicity,  only  turbulent  velocity  fluctuations  generated  by  ground-based 
wind  shear  are  explicitly  included  in  the  calculations;  thermal  turbulence  and  tur¬ 
bulence  generated  through  interactions  of  the  atmospheric  flow  with  the  barrier  are 
neglected.  This  idealization  is  reasonable  when  the  wind  blows  parallel  to  the  barrier. 
For  oblique  wind  angles,  one  can  expect  increased  turbulence  downwind  of  the 
barrier.  The  strength  of  the  turbulent  velocity  fluctuations  is  specified  with  the 
friction  velocity  w*,  which  is  related  to  the  dissipation  rate  through  the  equation 
£  rr  wj/0.4z,  where  z  is  set  to  2.8  m  [16].  (The  turbulence  parameters  Ly  and  e  could 
have  been  made  height  dependent  in  generating  the  QW  field,  but  this  would  com¬ 
plicate  interpretation  of  the  results.)  Furthermore,  only  scattering  elements  with  a 
line-of-sight  to  both  the  source  and  receiver  are  considered.  Because  the  scattering 
model  (based  on  von  Karman’s  spectrum)  is  valid  in  the  forward  directions,  we  need 
not  restrict  ourselves  to  scattering  angles  greater  than  20'’,  as  was  necessary  for 
earlier  calculations  [1]  based  on  Kolmogorov’s  spectrum.  The  reason  for  considering 
only  the  line-of-sight  paths  is  to  make  a  clear  physical  distinction  between  the  dif¬ 
fraction  process  and  scattering  by  intrinsic  atmospheric  turbulence.  Scattering  along 


Fig.  8.  Comparison  of  scattering  calculations  performed  with  the  Kolmogorov,  von  Karman,  and  QW 
models.  A  total  of  128  QW  size  classes,  ranging  from  ai  =  1  m  to  =  0.025  m,  were  used. 
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non-line-of-sight  paths  likely  involves  complicated  interactions  between  the  wind 
field  and  barrier,  which  have  the  additional  effect  of  refracting  sound.  The  reader  is 
referred  to  [22,23]  for  discussions  and  models  of  these  phenomena. 

The  calculated  spatial  extent  of  the  scattering  volume  for  this  geometry  at 
/  =  2000  Hz  is  shown  in  Fig.  6.  This  figure  was  created  by  calculating  the  ) 
corresponding  to  each  dVi  (using  Eq.  (43)),  sorting  the  results  in  order  of  increasing 
(\pff),  and  then  cumulatively  summing  the  resulting  sequence.  Note  that  the 
strongest  scattering  occurs  in  a  small  region  immediately  above  the  barrier.  Scat¬ 
tering  weakens  substantially  a  few  meters  above  the  barrier  or  5-10  m  toward  the 
source  or  receiver.  Results  for  other  frequencies  (not  shown)  demonstrate  that  the 
appearance  of  the  scattering  volume  does  not  change  substantially  over  the  fre¬ 
quency  range  considered  in  this  paper. 

Sound  pressure  levels  (SPL)  of  the  scattered  fields  corresponding  to  the  von 
Karman  (Eq.  (43)),  Kolmogorov  (large  kLy  limit  of  Eq.  (43)),  and  QW  models  are 
plotted  in  Fig.  7.  The  SPL  is  represented  here  as  101og(|psP)/bof,  where  po  is  the 
pressure  field  in  free  space.  For  comparison  purposes,  a  similar  plot  is  presented  in 
Fig.  8  except  with  the  size  of  the  largest  QW  class  reduced  fivefold  and  that  of  the 
smallest  QW  class  increased  fivefold  (ai  =  1  m  and  =  0.025  m,  with  Ly  in  the  von 
Karman  model  changed  accordingly).  At  low  frequencies,  the  von  Karman  and  QW 
models  give  essentially  the  same  result.  This  is  to  be  expected,  since  they  produce  the 


Fig.  9.  Contribution  to  the  scattered  mean-square  pressure  as  a  function  of  eddy  size, 
Results  for  QW  scattering  at  the  frequencies  500,  2000,  and  8000  Hz  are  shown. 
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same  turbulence  spectrum  within  the  energy-containing  and  inertial  subranges  (Fig. 
2).  The  Kolmogorov  spectrum  predicts  a  higher  scattered  SPL  at  low  frequency, 
because  it  lacks  a  realistic  spectral  roll-off  at  low  wavenumber.  This  effect  is  barely 
evident  in  Fig.  7  but  is  quite  clear  in  Fig.  8,  due  to  the  smaller  value  of  a|.  At  higher 
frequency,  the  Kolmogorov  and  von  Karman  models  give  the  same  results,  whereas 
the  QW  model  predicts  a  somewhat  lower  SPL.  This  is  a  result  of  the  spectral  roll-off 
in  the  dissipation  subrange,  which  is  included  in  the  QW  model  but  not  the  others. 
One  may  deduce  from  Fig.  7  that  the  Kolmogorov  spectrum  is  a  reasonable  ap¬ 
proximation  for  the  scenario,  even  at  frequencies  as  low  as  250  Hz.  This  outcome  is 
very  dependent  upon  the  barrier  geometry  and  outer  length  scale  of  the  turbulence, 
however. 

Fig.  9  shows  the  scale  of  the  eddies  (QWs)  responsible  for  the  scattering  for  the 
frequencies  500,  2000,  and  8000  Hz.  Specifically,  the  ratio  {bsP)/(IPsP)  at  the  re¬ 
ceiver  position  is  plotted.  The  maxima  for  these  curves  appear  at  —  0.30  m  for  500 
Hz,  fla  0.080  m  for  2000  Hz,  and  0.022  m  for  8000  Hz.  As  shown  in  Section  3, 
the  scattering  cross-section  is  maximized  for  the  size  class  a  such  that 
aala\  =  y/\l / 6/[2A:Lv  sin (0/2)] .  Taking  3  m  as  the  nominal  height  of  the  scattering 
volume,  and  assuming  that  the  horizontal  center  of  the  scattering  volume  is  above 
the  barrier,  implies  a  scattering  angle  of  about  20°  for  the  unreflected  path  SOR. 
With  ax  we  therefore  have  c^260/’'  m  for  SOR.  Maximum  scattering  is 


Fig.  10.  Scattered,  diffracted,  and  total  SPL  as  a  function  of  frequency.  The  barrier  height  is  2.5  m  and  the 
friction  velocity  (a  measure  of  the  strength  of  the  turbulence)  is  0.6  m  s"'. 
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thus  expected  to  occur  at  «« 0.52  m  for  500  Hz,  0.13  m  for  2000  Hz,  and 
aoc  ^  0.032  m  for  8000  Hz.  These  values  are  45-70%  larger  than  the  locations  of  the 
calculated  peaks.  However,  it  must  be  kept  in  mind  that  the  the  actual  scattering 
volume  is  distributed  in  space,  and  that  the  regions  distant  from  the  barrier  are 
characterized  by  higher  scattering  angles.  Furthermore,  the  ground-reflected  paths 
have  higher  scattering  angles.  A  repeat  of  this  exercise  for  the  double  reflection  path, 
WOR,  yields  a  scattering  angle  of  about  35"^,  ~  150/'“*  m,  Oa  0.30  m  for  500  Hz, 

fla  0.075  m  for  2000  Hz,  and  ««  0.019  m  for  8000  Hz.  These  values  are  much 

closer  to  the  theoretical  result,  thereby  suggesting  the  correctness  of  the  theory.  It  is 
also  interesting  that  the  scattering  contributions  have  a  distinctly  normal  (bell-curve) 
appearance  when  plotted  as  a  function  of  the  logarithm  of  the  eddy  size.  Most  of  the 
scattering  occurs  from  eddy  sizes  between  one-half  and  twice  the  peak  scattering  size. 

The  frequency  dependence  of  the  scattered,  diffracted,  and  total  received  SPL  is 
shown  in  Fig.  10.  The  diffracted  field  was  calculated  with  Pierce’s  formula  [20,  Eq. 
(9-9.1)].  The  scattered  field  in  this  example  was  calculated  from  the  friction  velocity 
w*  ==  0.6  m  s“*,  which  corresponds  to  Cl  =  0.47  m'^/^s"^  at  2.8  m  height  and  is  a 
rather  large  value  for  the  atmosphere  [24].  (In  practice,  the  barrier  may  generate 
additional  turbulence,  effectively  raising  m*.)  Only  for  frequencies  above  approxi¬ 
mately  4  kHz  does  scattering  significantly  affect  the  observed  SPL.  This  outcome  is  in 


Fig.  11.  Scattered,  diffracted,  and  total  SPL  as  a  function  of  barrier  height.  The  heights  and  horizontal 
distances  of  the  source  and  receiver  are  held  constant.  The  calculation  is  for  /  =  8000  Hz  and  a  friction 
velocity  of  0.6  m  s”L 
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broad  agreement  with  Fig.  8  from  [1],  although  the  details  of  the  frequency  depen¬ 
dence  are  different  here  due  to  the  inclusion  of  the  ground-reflection  paths. 

Fig.  1 1  shows  the  received  SPL  at  8000  Hz  as  a  function  of  the  barrier  height.  The 
heights  of  the  source  and  receiver,  and  their  horizontal  distances  from  the  barrier, 
were  held  fixed  to  the  values  given  earlier.  For  barrier  heights  less  than  about  4  m, 
the  scattered  contribution  is  roughly  equal  to  the  diffracted  one.  At  greater  heights, 
and  therefore  at  larger  scattering  angles,  the  diffracted  contribution  is  more  prom¬ 
inent.  Fig.  12  is  similar  to  Fig.  11,  except  that  the  received  SPL  at  8000  Hz  is  shown 
as  a  function  of  the  horizontal  separation  between  the  barrier  and  receiver  (jCr).  The 
barrier  height  for  Fig.  12  is  again  2,5  m  and  the  horizontal  distance  to  the  source 
remains  fixed  at  Xs  =  — 10  m.  As  the  distance  of  the  receiver  from  the  barrier  in¬ 
creases,  the  scattering  angle  becomes  smaller  and  therefore  the  scattered  energy 
becomes  more  significant.  Also  note  that  as  the  distance  increases  the  effects  of 
multiple  scattering  in  the  direction  of  sound  propagation  become  more  important. 
These  effects,  however,  result  mainly  in  fluctuations  in  the  phase  and  amplitude  of  a 
sound  wave  reaching  the  receiver  while  the  intensity  of  the  wave  remains  approxi¬ 
mately  unchanged. 

Last,  Fig,  13  shows  the  scattered,  diffracted,  and  total  received  SPL  at  8000  Hz  as 
a  function  of  the  friction  velocity.  As  increases  past  0.6  m  s“',  the  scattered 


Fig.  12.  Scattered,  diffracted,  and  total  SPL  as  a  function  of  horizontal  distance  between  the  barrier  and 
receiver.  The  heights  of  the  source  and  receiver,  and  the  horizontal  distance  to  the  source,  are  held  con¬ 
stant.  The  calculation  is  for  /  =  8000  Hz,  a  barrier  height  of  2.5  m  and  a  friction  velocity  of  0.6  m  s“'. 
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Fig.  13.  Scattered,  diffracted,  and  total  SPL  as  a  function  of  the  friction  velocity  of  the  turbulence.  Results 
for  the  friction  velocity  varying  from  0.1  to  1  m  s“*  are  shown.  The  calculation  is  for  f  =  8000  Hz  and  a 
barrier  height  of  2.5  m. 


contribution  overtakes  the  diffracted  one.  Taken  together,  Figs.  10-13  demonstrate 
that  turbulent  scattering  is  significant  primarily  at  high  frequencies  and  for  relatively 
short  barriers.  This  outcome  is  consistent  with  Fig.  4,  which  shows  that  the 
scattering  cross-section  is  highest  for  nearly  forward  scattering  at  large  values  of  kLy, 
Strong  turbulence,  either  already  intrinsic  to  the  atmosphere  or  injected  by  flow 
disturbances  from  the  barrier  or  other  nearby  objects,  must  be  present  for  significant 
scattering. 


5.  Conclusion 

As  was  first  demonstrated  by  Daigle  [1],  turbulence  diminishes  the  effectiveness  of 
noise  barriers  by  scattering  sound  energy  over  them.  Daigle’s  calculation  of  the 
scattered  sound  energy  was  based  on  the  scattering  cross-section  for  inertial  sub¬ 
range  turbulence,  namely  the  large  kLy  limit  of  Eq.  (44).  In  this  paper,  we  have 
performed  scattering  calculations  similar  to  Daigle’s,  although  the  recently  intro¬ 
duced  quasi-wavelet  (QW)  representation  was  used.  The  QW  representation  mimics 
the  conceptual  picture  of  turbulence  consisting  of  a  number  of  discrete,  randomly 
placed  eddies  whose  sizes  span  a  range  of  spatial  scales.  Parameters  of  the  QW 
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representation,  such  as  the  number  of  eddies,  their  size  distribution,  and  angular 
velocities,  can  be  chosen  to  reproduce  known  models  for  the  atmospheric  turbulence 
spectrum.  When  applied  to  scattering,  the  QW  representation  provides  a  concise 
decomposition  of  the  scattering  as  a  function  of  the  eddy  size  and  location.  Also, 
since  the  QW  parent  function  was  derived  in  this  paper  on  the  basis  of  von  Karman’s 
spectrum,  it  can  be  applied  to  scattering  from  the  energy-containing  subrange  and 
therefore  to  small  scattering  angles.  ^ 

An  expression  for  the  scattering  cross-section  of  an  ensemble  of  randomly  ori¬ 
ented  QWs  of  a  particular  size  was  derived  in  this  paper.  When  the  scattered  energies 
of  all  of  the  QW  size  classes  were  added  together,  and  the  continuous  size  limit  taken, 
the  total  scattering  cross-section  agreed  with  a  previously  derived  result  for  the  von 
Karman  turbulence  spectrum.  An  expression  for  the  eddy  size  that  maximizes 
scattering,  as  a  function  of  the  acoustic  wavenumber  and  scattering  angle,  was  also 
derived.  Calculations  showed  that,  for  a  typical  barrier  condition,  most  of  the 
scattered  energy  originates  from  eddies  in  the  size  range  of  approximately  one-half  to 
twice  the  size  of  the  eddies  responsible  for  maximum  scattering. 

The  results  in  this  paper  suggest  that  scattering  behind  the  barrier  attributable  to 
“line-of-sight”  eddies,  namely  those  with  a  line  of  sight  to  both  the  source  and  re¬ 
ceiver,  is  generally  significant  only  for  frequencies  above  several  kilohertz,  for 
sources  and  receivers  no  more  than  a  few  meters  below  the  top  of  the  barrier,  and  for 
very  turbulent  atmospheric  conditions. 
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